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Constraints on dark energy from current observational data are sensitive to how distances are 
measured from Type la supernova (SN la) data. We find that flux- averaging of SNe la can be used 
to test the presence of unknown systematic uncertainties, and yield more robust distance measure- 
ments from SNe la. We have applied this approach to the nearby+SDSS+ESSENCE+SNLS+HST 
set of 288 SNe la, and the "Constitution" of set 397 SNe la. Combining the SN la data with cos- 
mic microwave background anisotropy data from Wilkinson Microwave Anisotropy Probe five year 
observations, the Sloan Digital Sky Survey baryon acoustic oscillation measurements, the data of 
69 gammay-ray bursts, and the Hubble constant measurement from the Hubble Space Telescope 
project SHOES, we measure the dark energy density function X(z) = px(z)/px(0) as a free func- 
tion of redshift (assumed to be a constant at z > 1 or z > 1.5). Without flux-averaging of SNe 
la, the combined data using the "Constitution" set of SNe la seem to indicate a deviation from a 
cosmological constant at ~95% confidence level at < z < 0.8; they are consistent with a cosmo- 
logical constant at ~68% confidence level when SNe la are flux-averaged. The combined data using 
the nearby+SDSS+ESSENCE+SNLS+HST data set of SNe la are consistent with a cosmological 
constant at 68% confidence level with or without flux-averaging of SNe la, and give dark energy 
constraints that are significantly more stringent than that using the "Constitution" set of SNe la. 
Assuming a flat universe, dark energy is detected at > 98% confidence level for z < 0.75 using the 
combined data with 288 SNe la from nearby+SDSS+ESSENCE+SNLS+HST, independent of the 
assumptions about X(z > 1). We quantify dark energy constraints without assuming a flat universe 
using the dark energy Figure-of-Merit (FoM) for both X(z) and a dark energy equation-of-state 
linear in the cosmic scale factor. 
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I. INTRODUCTION 



The observational evidence of cosmic acceleration [J, 2| 
continues to strengthen with time (e.g., [E 0, [E [3 0, 
Sllll)- However, we are still in the dark about the 
nature of the observed cosmic acceleration; whether it is 
due to an unknown energy component [ll|, EE EE EE 
EE EE EE Ell, i- e -, dark energy, or a modification of 



general relativity [1 
recent reviews, see [2 
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123, I23|, J2J 



For 



Differentiating dark energy from modified gravity as 
the cause for the observed cosmic acceleration requires 
ambitious observational projects of galaxy redshift sur- 
veys |33| and weak lensing surveys [341 ] that are still in 
the planning stages. Currently available data allow us to 
begin to test whether the observed cosmic acceleration 
could be due to a cosmological constant, the simplest 
possibility for dark energy. 

A critical challenge in our understanding of dark en- 
ergy is the control of known and unknown systematic 
uncertainties. Here we investigate methods for measur- 
ing distances from Type la supernova (SN la) data, and 
their impact on dark energy constraints from current ob- 
servational data. 

In the usual likelihood analysis of distances to SNe la 
that directly utilizes the measured brightness of each SN 
la, the statistics is dominated by the SNe la with the 



smallest measurement errors. This could lead to biased 
distance measurements in the presence of unknown sys- 
tematic errors. We will show that this pitfall can be 
removed by flux-averaging SN la data in redshift bins 
(which also reduces the bias in estimated parameters due 
to weak gravitational lensing of the highest redshift SNe 
la). 

We describe our method in Sec. II, present our results 
in Sec. Ill, and conclude in Sec. IV. 



II. METHOD 

The minimal way to include dark energy in the stan- 
dard cosmological framework is to add a new energy com- 
ponent with density px(z). The Friedman equation be- 
comes 



= nl [fi ro (i + zf + n r (i + zf + n fe (i + zf + n x x{ z 

where fl m + fi r + Qk + &x — 1, and the dark energy 
density function X(z) is defined as 



X(z) 



Px{z) 



(2) 
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Note that fl r <C f2 m , thus the fl r term is usually omitted 
in dark energy studies, since dark energy should only be 
important at late times. 

The comoving distance to an object at redshift z is 
given by: 

r(z) =cff - 1 Ififcl-^sinnllfifcl 1 / 2 T(z)}, (3) 
r(z) = r^)' E{z)=H{z)/H <i 

where sinn(a;) = sin(x), x, sinh(x) for 0,^ < 0, Ct^ = 0, 
and > respectively. 



A. Analysis of SN la Data 

The published SN la data sets usually give the distance 
modulus to each SN la (marginalized over calibration 
parameters): 



fiQ — m — M = 5 log 



Mpc 



25, 



(4) 



with 



where the luminosity distance d^{z) = (1 + z) r(z), 
the comoving distance r(z) given by Eq.Q. 

We consider the two latest compilations of 
SN la data, the "Constitution" set of 397 
SNe la by Hicken et al. (2009) [HI, and the 
nearby+SDSS+ESSENCE+SNLS+HST set of 288 
SNe la by Kessler et al. (2009) [36|. The "Constitution" 
set used SALT (37[ to fit the SN la light curves. We use 
the nearby+SDSS+ESSENCE+SNLS+HST set that 
used SALT2 [H for SN la lightcurve-fitting. We do not 
use the nearby+SDSS+ESSENCE+SNLS+HST set that 
used MLCS [ai HI Ell for SN Ia lightcurve-fitting, as 
the MLCS method appears to introduce some systematic 
biases [36j . 

It has been noted that the "Constitution" set of SNe 
Ia (and other prior compilations of SNe Ia) appear to 
be inhomogeneous, possibly due to unknown systematic 
effects [42|, |43j . Here we use the flux-averaging of SNe Ia 
to help reduce the impact of unknown systematic effects. 

Flux-averaging of SNe Ia was proposed to reduce the 
effect of the weak lensing of SNe Ia on cosmological pa- 
rameter estimation The basic idea is that because 
of flux conservation in gravitational lensing, the average 
magnification of a large number of SNe Ia at the same 
redshift should be unity. Thus averaging the observed 
flux from a large number of SNe Ia at the same redshift 
can recover the unlensed brightness of the SNe Ia at that 
redshift. 

Wang & Mukherjee (2004) [45] and Wang (2005) 
[HI developed a consistent framework for flux-averaging 
SNe Ia. Appendix A of Wang & Mukherjee (2007) 
[47l | describes in detail the recipe for flux-averaging 
SNe Ia used here; a public code i s available at 
|http://www.nhn.ou.edu/~wang/SNcode/[ Here we do 
not marginalized over i?o, since it is absorbed in the "nui- 
sance parameter" Mgjv that is marginalized over in our 
likelihood analysis. 



Since the SNe Ia in each redshift bin are flux-averaged 
and not used directly in the likelihood analysis, the sys- 
tematic bias that results from unknown systematic errors 
from individual SNe Ia can be minimized. For homo- 
geneous data without unknown systematic errors, flux- 
averaging should not change the results qualitatively at 
intermediate redshifts (where the lensing effect is negli- 
gible). Thus in addition to reducing lensing and lensing- 
like systematic effects, flux-averaging of SN Ia data pro- 
vide a useful test for the presence of unknown system- 
atic effects at intermediate redshifts. We limit the flux- 
averaging to SNe Ia at z > 0.2, since the SNe Ia at lower 
redshifts are better understood. 



B. CMB data 

CMB data give us the comoving distance to the 
photon-decoupling surface r(z*), and the comoving sound 
horizon at photo-decoupling epoch r s (z») [48]. Wang & 
Mukherjee 2007 H3 showed that the CMB shift param- 
eters 



R 



)Jtt m H$r(z*), l a =7rr(z*)/r a (z*), 



(5) 



together with ujf, = Qbh 2 , provide an efficient summary 
of CMB data as far as dark energy constraints go. Us- 
ing {R 7 l a ,uJb} is equivalent to using {R,l a ,z*} as CMB 
distance priors [Hoj ]. 

The comoving sound horizon at redshift z is given by 



r s (z) 



* c, dt' 



cH n 



dz 



E(z>y 



da' 



3(1 + # 6 a')a' 4 E 2 {z') 



(6) 



where a is the cosmic scale factor, a = 1/(1 + z), 
and a 4 E 2 (z) = Q m (a + a oq ) + Sl^a 2 + QxX(z)a 4 , 
with a oq = firad/^m = 1/(1 + Zeq), and z cq = 
2.5 x lO i n m h 2 (T C MB/2.7K)- i . The sound speed is 

c s = 1/^3(1 + R b a), with R b a = 3pb/(4p 7 ), Rb = 

31500fW(T CM B/2.7Kr 4 . We take T C mb = 2.725 fol- 
lowing Komatsu et al. (2009) since we will use the 
CMB bounds derived by them. 

Here we use the covariance matrix of 
\R(z*), lg,(z*), z*,r 8 (zd)] from the five year WMAP 
data [3j, [49J , which includes the comoving sound horizon 
at the drag epoch r s (zd). Note that z* is given by the 
fitting formula (Hlj : 



= 1048 [1 + 0.00124(fM 2 )- - 738 ] [I + gi(Q. m h 



2\92] 



where 



91 

.92 



0.0783 (Q b h 



2\-0.238 



1 + 39.5 (fW) - 763 

0.560 
l + 21.1(fi 6 /i 2 ) 1 - 81 



(7) 

(8) 
(9) 
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The redshift of the drag epoch Zd is well approximated 

by m 



Zd = 



1291(r2 m /i 



2\0.251 



1 + 0.659(f7 m /i 2 ) 0828 



2\fc21 



(10) 



where 



h = 0M3(n m h 2 y 0A19 [l + 0.607(O„/i 2 ) - 674 ] ,(11) 
b 2 = 0.238(a«/i 2 ) - 223 - (12) 

CMB data are included in our analysis by adding the 
following term to the y 2 of a given model with p\ — 

R(z*), p 2 = l a (z*), P3 = and p A = r s {z d ): 

Xcmb = &Pi [Cov cM B (PhPj)] A Pn &Pi = P*-P? ata , 

(13) 

where pd-ata are ^ ne maximum likelyhood val- 
ues, and Covcmb is the covariance matrix of 
[R{z*),l a (z*),z*,r s {z d )] BIH. 



where r(z) is the comoving distance at z. 

The GRB data are included in our analysis by adding 
the following term to the \ 2 of a given model: 



Xgrb = [ A r p (zi)} ■ (Cov^g) ■ [Ar p (zj)] 



Ar^(zi) = r^ iata (zi)-r^(zi) 



(17) 



where r p (z) is defined by Eq. (|16p . The covariance matrix 
is given by 



(Covg^b),., = v{rp{zi))o(rp(zj)) (Cov 



GRB) 



(18) 



where Govgrb is the normalized covariance matrix from 
Table 3 of Wang (2008b) 0, and 

<rfc(Zi)) = <j(r^(z l )) + , if r-(z) > ~ v (z) data ; 
a{r^{ Zi )) = a(r^(zi))-, if rj(z) < r>(z) data ,(19) 

where a (r^(zi)) + and a (fp~(zi))~ are the 68% C.L. errors 
given in Table 2 of Wang (2008b) [H. 



C. Baryon Acoustic Oscillation Data 



E. Dark energy parametrization 



For the baryon acoustic oscillation (BAO) data, we use 
the measurement of d z = r s (z d )/ Dy(z) at z = 0.2 and 
z = 0.35 by Pcrcival ct al. (2009) [11], where 



D V (z) EE 



"(z) 2 cz 



1/3 



(14) 



The inverse covariance matrix of (do. 2, do. 35), Cov s ^ , 
is given by [53[: Cov B ^ n = 30124, Cov^ 012 
1 7227. and Cov B \ 22 = 86977. 

BAO data are included in our analysis by adding the 
following term to the \ 2 of a given model with p\ = di 
and pi = do.35: 



0.2 



XbAO = A Pi [ Gov BAo(.PiiPj)] A Pj> A P* =Pi~P' 



data 



where pf a ta 
0.109711. 



(15) 

l$ a 2 ta = 0.1905, and pf ta = d$ a $ = 



D. Gammay-ray Burst Data 

We add gammay-ray burst (GRB) data to our analy- 
sis, since these are complementary in redshift range to the 
SN la data. We use GRB data in the form of the model- 
independent GRB distance measurements from Wang 
(2008b) [HI, which were derived from the data of 69 
GRBs with 0.17 < z < 6.6 from Schaefer (2007) [H|. 

The GRB distance measurements are given in terms of 

M 



~( z i) 



r p(z) 
r p (0.17) : 



(1 + zfl 2 H 



ch 



r(z), 



(16) 



Since we are ignorant of the the true nature of dark en- 
ergy, it is useful to measure the dark energy density func- 
tion X(z) = px(z)/ px(0) as a free function of redshift 
[5(|[57], 58]. Here we parametrize X(z) by cubic-splining 
its values at z = 0.25, 0.5, 0.75, and 1.0, and assume that 
X(z > 1) = X(z = 1). Fixing X(z > 1) reflects the limit 
of current data, and avoids making assumptions about 
early dark energy that can be propagated into artificial 
constraints on dark energy at low z [4?], H3] ■ 

For comparison with the work of others, we also con- 
sider a dark energy equation of state linear in the cosmic 
scale factor a, iux(a) = wq + (1 — a)w a [5S^ . A related 
parametrization is [501 ] 



w x (a) 



Or 



a c - 
a c w 



w 



W0.5 



1 



a(w . 5 



w .5 



Wq) 



1 



(20) 



with a c = 2/3 (i.e., z c = 0.5), and W0.5 — wx(z = 0.5). 
Eq. (|2"0"|) is related to wx(z) — wq + (1 — a)w a by setting 



Wq.5 ~ wo 



1 



or w . 5 = wq + (1 - a c )w a - 

(21) 

Wang (2008a) [50j showed that (wo, W0.5) are much less 
correlated than (wq, w a ), thus are a better set of param- 
eters to use. 



III. RESULTS 

We perform a Markov Chain Monte Carlo (MCMC) 
likelihood analysis 60] to obtain O(10 6 ) samples for 
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each set of results presented in this paper. We as- 
sum flat priors for all the parameters, and allow ranges 
of the parameters wide enough such that further in- 
creasing the allowed ranges has no impact on the re- 
sults. The chains typically have worst e-values (the vari- 
ance(mean)/mean(variance) of 1/2 chains) much smaller 
than 0.005, indicating convergence. The chains are sub- 
sequently appropriately thinned to ensure independent 
samples. 

In addition to the SN la, CMB, BAO, and GRB data 
discussed in SecllJ we impose a prior of H Q — 74.2 ± 
3.6kms -1 Mpc -1 , from the Supernovae and Hq for the 
Equation of State (SHOES) program on the HST [H|. 

We do not assume a flat universe unless specifically 
noted. In addition to the dark energy parameters de- 
scribed in Sec lIIEi we also constrain cosmological pa- 
rameters (fi m , fife, h, ujb, Msn), where uji, = flbh 2 , and 
Msn is used to model the absolute distance scale of the 
SNe la. 



Model-Independent Distance Measurements 
from SNe la 



For an intuitive understanding of the SN la data, we 
parametrize a scaled distance 



/ \ t(z) , , r 

cH n z 



(22) 



by its values at Zi = 0.125i, i — 1,2,. ..,10, and z\\ = 
1.551 (the highest redshift of SNe la in the "Constitu- 
tion" and nearby+SDSS+ESSENCE+SNLS+HST data 
sets), and measure r p (zi) (i — 1,2, ...,11) from the SN 
la data. The r p (z) at arbitrary z is given by cubic 
spline interpolation, thus no assumptions are made about 
cosmological models. The {r p (zi)} thus provide model- 
independent distance measurements from SNe la. 

Note that the scaled distances defined in Eqs. ([22|) and 
(fTS)) are similar, but with different choices of the power 
of {l + z) used; this choice is made to make the scaled dis- 
tance as flat as possible over the redshift range of interest, 
in order to maximize the accuracy for cubic-spline. 

For cosmological models allowed by current data, 
the accuracy of interpolating r p (z) using r p (zi) (i — 
0, 1, 2, 11, with z Q = and r p (0) = 1) is around 0.4%, 
compared to the exact distances. 

FigsfTJU show the distances mea- 
sured from the "Constitution" and 
nearby+SDSS+ESSENCE+SNLS+HST data sets. 
Clearly, the nearby+SDSS+ESSENCE+SNLS+HST 
set of 288 SNe la gives tighter constraints on dis- 
tances, although it contains 109 fewer SNe la than the 
"Constitution" set. 

We have flux-averaged the SNe la at z > 0.2, to 
minimize any likely systematic biases due to lensing 
or unknown systematic effects. It is interesting to 
note that flux-averaging has a more significant effect 
on the "Constitution" set. This is as expected, since 



n =.27, n„=.73 



w x =-0.9 
• w-=-l.l 



r p (z) . . I 




397 SNe la ("Constitution" set) 
\ SNe la flux-averaged with dz=0.03 for z£0.2 
© no flux-averaging 




FIG. 1: Distance measurements from the "Constitution" set of 397 
SNe la. and nearby+SDSS+ESSENCE+SNLS+HST data sets. 



n =.27, n„=.73 



w x =-0.9 
w,= -l.l 




288 SNe la (inc. SDSS) 

\ SNe la flux-averaged with dz=0.03 for z^0.2 
© no flux-averaging 




FIG. 2: Distance measurements from the 

nearby+SDSS+ESSENCE+SNLS+HST set of 288 SNe la. 



the "Constitution" set is less homogeneous than the 
nearby+SDSS+ESSENCE+SNLS+HST data set, which 
include 103 SNe la from SDSS at 0.04 < z < 0.42 
Flux-averaging brings both data sets closer to the pre- 
diction of a flat universe dominated by a cosmological 
constant. 



5 




FIG. 3: Dark energy density function X(z) = px (z) / px (0) mea- 
sured from combining SN la data (the "Constitution" set of 397 
SNe la) with CMB, BAO, GRB data, and imposing the SHOES 
prior on Hq. The 68% (shaded) and 95% confidence level regions 
are shown. 

Table I lists the distances measured from the 
nearby+SDSS+ESSENCE+SNLS+HST data set of 288 
SNe la, flux-averaged with dz = 0.03 at z > 0.2. Table 
II gives the corresponding normalized covariance matrix. 
Note that the distances measured from SNe la, r p (zi) 
(i = 1,2, ...,11), are only moderately correlated. Since 
the measurements of r p (Zi) (i — 1,2, ...,11), as given by 
Tables I and II, encode all the distance information from 
the nearby+SDSS+ESSENCE+SNLS+HST data set of 
288 SNe la independent of a cosmological model, they 
can be used in place of the full data set following the 
prescription given in Sec lIIDl 

B. Measurements of the Dark Energy Density 
Function and H(z) 

The most transparent way to see how current data 
compare to the prediction of the cosmological con- 
stant model is to measure the dark energy density 
function. FigOU] show the dark energy density func- 
tion measured from combining SN la data with CMB, 
BAO, GRB data, and imposing the SHOES prior 
on Hq, for SNe la from the "Constitution" and the 
nearby+SDSS+ESSENCE+SNLS+HST data sets re- 
spectively, without assuming a flat universe. Fig|5] and 
FigJS]are similar to Fig[3]and FigQJ except a flat universe 
is assumed. 

Again, the same trend as in FigsQ][2] is seen: using 
the nearby+SDSS+ESSENCE+SNLS+HST data set of 
SNe la gives much more stringent constraints on dark 
energy than using the "Constitution" set of SNe la, and 



FIG. 4: Dark energy density function X(z) = 

px( z )/px(0) measured from combining SN la data 
(nearby+SDSS+ESSENCE+SNLS+HST data set of 288 SNe la) 
with CMB, BAO, and GRB data, and imposing the SHOES prior 
on Hq. The 68% (shaded) and 95% confidence level regions are 
shown. 

gives measurements that are closer to a cosmological con- 
stant. Flux- averaging again has larger impact on the re- 
sults from using the "Constitution" set of SNe la, and 
brings the measurements closer to that predicted by a 
cosmological constant. 

Note that the X(z) measured using the "Constitution" 
set of SNe la actually deviates from a cosmological con- 
stant at ~2<t without flux-averaging; the downturn in the 
measured X(z) at z ~ 0.25 — 0.5 is consistent with that 
found by Huang et al. (2009) [g^]. Since this apparent 
variation in X{z) disappears in the results from using 
the nearby+SDSS+ESSENCE+SNLS+HST set of SNe 
la without flux-averaging, it is likely that it originated 
from unknown systematic effects in the "Constitution" 
set of SNe la. 

FiglSJa) shows that assuming a flat universe, dark 
energy is detected at > 2a at z = 1 using the 
nearby+SDSS+ESSENCE+SNLS+HST data set of 288 
SNe la, together with CMB, BAO, and GRB data. 
FigEfb) is the same as Fig^a), but with X(z = 1.5) 
added as a parameter, and assuming X(z > 1.5) = 
X(z = 1.5). The X(z) results at z < 0.75 remain about 
the same, while the errors for X(z — 1) are larger (but 
the mean value oiX(z) remains about the same) because 
of the additional parameter allowed. Independent of the 
assumptions about X[z) at z > 1, dark energy is detected 
at > 98% confidence level for z < 0.75. 

Figs OS] are the measurements of the cosmic expansion 
history H(z) that correspond to Fig|3][4] 

Table III gives the dark energy density func- 
tion X(z) = px{z)l px(0) and the cosmic ex- 
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TABLE I: Distances measured from 288 SNe la with 68% C.L. upper and lower uncertainties. 





z 


TZ data (z) 
p \ ) 


\ V) 


<j(rJz)) + 
V P \ } t 


V P\ ) ) 


o 


0.0 


1.0000 








1 


0.125 


1.0044 


0.0129 


0.9919 


1.0171 


2 


0.25 


1.0523 


0.0129 


1.0396 


1.0651 


3 


0.375 


1.0292 


0.0199 


1.0095 


1.0490 


4 


0.5 


1.0638 


0.0244 


1.0397 


1.0878 


5 


0.625 


1.0311 


0.0245 


1.0067 


1.0556 


6 


0.75 


1.0895 


0.0321 


1.0575 


1.1212 


7 


0.875 


1.0160 


0.0356 


0.9807 


1.0512 


8 


1.0 


1.0629 


0.0398 


1.0234 


1.1022 


9 


1.125 


1.0725 


0.0747 


0.9996 


1.1452 


10 


1.25 


1.0557 


0.0427 


1.0136 


1.0977 


11 


1.551 


0.9682 


0.1505 


0.8203 


1.1151 



TABLE II: Normalized covariance matrix of distances measured from 288 SNe la 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


1 


1.0000 


0.2539 


0.5219 


0.2016 


0.3354 


0.2148 


0.2030 


0.1741 


0.0873 


0.1728 


0.0349 


2 


0.2539 


1.0000 


0.1386 


0.2120 


0.1496 


0.1483 


0.1068 


0.1192 


0.0635 


0.0971 


0.0323 


3 


0.5219 


0.1386 


1.0000 


0.1288 


0.2575 


0.1413 


0.1421 


0.1130 


0.0656 


0.1106 


0.0350 


4 


0.2016 


0.2120 


0.1288 


1.0000 


-0.0694 


0.1473 


0.0235 


0.0917 


0.0280 


0.0650 


0.0117 


5 


0.3354 


0.1496 


0.2575 


-0.0694 


1.0000 


0.1236 


0.1273 


0.0677 


0.0537 


0.0809 


0.0234 


6 


0.2148 


0.1483 


0.1413 


0.1473 


0.1236 


1.0000 


-0.1503 


0.1709 


-0.0173 


0.0521 


-0.0163 


7 


0.2030 


0.1068 


0.1421 


0.0235 


0.1273 


-0.1503 


1.0000 


-0.0877 


0.0796 


0.0609 


0.0515 


8 


0.1741 


0.1192 


0.1130 


0.0917 


0.0677 


0.1709 


-0.0877 


1.0000 


0.0811 


0.0299 


-0.0129 


9 


0.0873 


0.0635 


0.0656 


0.0280 


0.0537 


-0.0173 


0.0796 


0.0811 


1.0000 


0.1779 


0.3163 


10 


0.1728 


0.0971 


0.1106 


0.0650 


0.0809 


0.0521 


0.0609 


0.0299 


0.1779 


1.0000 


-0.4585 


11 


0.0349 


0.0323 


0.0350 


0.0117 


0.0234 


-0.0163 


0.0515 


-0.0129 


0.3163 


-0.4585 


1.0000 



pansion history H(z) measured from current data 
(nearby+SDSS+ESSENCE+SNLS+HST data set of 288 
SNe la, together with CMB, BAO, GRB data, and im- 
posing the SHOES prior on Ho). The H(z) measure- 
ments are derived using Eq.{T]). Tables IV and V give 
the normalized covariance matrices of the X(z) and H (z) 
measurements. Note that both the X(z) and H(z) mea- 
surements are only weakly correlated. 



C. Constraints on (wo,w a ), (u>o,uio.s) 

For comparison with the work of others, FigslMTUl 
show the 68% and 95% joint confidence level contours 
of (wo,w a ) and (100,11)0.5), from combining SN la data 
with CMB, BAO, GRB data, and imposing the SHOES 
prior on Ho, for SNe la from the "Constitution" and 
the nearby+SDSS+ESSENCE+SNLS+HST data sets 
respectively, 

It is interesting to note that flux-averaging has 



a larger impact on the (wo,w a ) and (wo^o.s) er- 
ror contours from the combined data using the 
nearby+SDSS+ESSENCE+SNLS+HST set of SNe la, 
compared to the combined data using the "Constitu- 
tion" of SNe la. This is in contrast to what we found 
when measuring the dark energy density function X(z) 
at z = 0.25,0.5,0.75,1 (see FigsOE]). A comparison of 
the correlation coefficients in Tables IV and VI reveal 
that assuming a linear dark energy equation-of-state re- 
sults in much stronger correlations among the dark en- 
ergy parameters. This obsures some of the information 
about dark energy contained in the data. 
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TABLE III: X(z), H(z), and cosmological parameters estimated from current data with 68% C.L. upper and lower uncertainties. 
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0.942 


0.059 
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X(z 


— 1 ft) 
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0.012 
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LOb 
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0.00062 
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0.02441 
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1.10 
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1.12 


H(z 


= 0.5) 


1.28 


0.043 


1.24 


1.32 


H(z 


= 0.75) 


1.51 
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1.37 
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1.66 


0.215 


1.45 
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TABLE IV: Normalized covariance matrix of X(z) from current data 
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2 
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1.0000 


-0.1273 


-0.0668 


0.0078 
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-0.1273 


1.0000 


-0.0374 


-0.2789 


3 


-0.0668 


-0.0374 


1.0000 


-0.3638 


4 


0.0078 


-0.2789 
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D. Dark Energy Figure of Merit from Current 
Data 

Wang (2008a) [5(| defined a general dark energy 
Figure-of-Merit (FoM) 

FoM = 1 (23) 

VdetCovai,^,.../*) 

where (/i, f%, .../jv) is the set of parameters that have 
been chosen to parametrize dark energy. 

The Dark Energy Task Force (DETF) defined the dark 
energy FoM to be the inverse of the area enclosed by 
the 95% confidence level contour of (wo,w a ) 63]. The 
areas enclosed by contours are difficult to calculate for 
real data, as these contours can be quite irregular (see 
FigsMTUl) . The definition of Eq. has the advantage of 
being easy to calculate for either real or simulated data. 
For (A,/ 2 ) = (w ,w a ), FoM of Eq. (f2"3")) is proportional 
to the FoM defined by the DETF for ideal Gaussian- 
distributed data, and the same as the relative FoM used 
by the DETF in Fisher matrix forecasts. 

Table VI shows the dark energy FoM from SN la data 
together with CMB, BAO, GRB data, and imposing the 
SHOES prior on for the dark energy density function 
X(z) = p x {z)/px(0) measured at z = 0.25,0.5,0.75,1.0 
(with X(z) given by cubic spline elsewhere and assum- 
ing X(z > 1) = X(l), see Sec lII E[) and a dark energy 



equation-of-state linear in a parametrized by (wq,wq^) 
and (wo,w a ) respectively. 

It is interesting to note that for {X(zi)} 
(zi=0.25, 0.5, 0.75, 1.0), the factor of improve- 
ment in the dark energy FoM is 5.2 when the 
nearby+SDSS+ESSENCE+SNLS+HST data set of 
288 SNe la is used instead of the "Constitution" set of 
397 SNe la (see Table VI); this quantifies the difference 
between Fig0] and Fig[3J For a linear dark energy 
equation of state, the improvement in FoM is a little 
less than ^50% (see Table VI). The same equation of 
state linear in a can be written in either (wq,wq^) or 
(wQ,w a ) (see Sec lII E| . As expected, the parameters 
(wq,wo.5) are significantly less correlated than (wo,w a ) 
|5fj |. thus represent a better choice of parameters. The 
improvement in the dark energy FoM is slightly larger 
for (wo, 100.5) than for (wo,w a ). 



E. The Impact of Adding CMB Constraints on 

r s (z d ) 

Compared to previous work using CMB distances pri- 
ors, we have added the constraints on r s {zd)- Not sur- 
prisingly, this has negligible impact on X(z) measured 
from current data, since we have assumed that X(z > 
l)=X(z=l). 

Figs fTTIIT2l show the difference in the marginalized 
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TABLE V: Normalized covariance matrix of H(z) from current data 





1 


2 


3 
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1.0000 


-0.1665 


-0.1212 
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1.0000 
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1.0000 
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TABLE VI: Dark energy FoM from current data 



data: CMB+BAO+GRB+ 


FoM f .({Z i }) 


a(w ) cr(wo.s) r WOlW05 FoMr-(iuo, wo.s) 


a(w ) a(w a ) r WOyWa FoM r (w ,w a ) 


288 SNe la (inc. SDSS) 
397 SNe la ("Constitution") 


230.7 
44.1 


0.1473 0.1688 -0.6676 54.0 
0.1841 0.2252 -0.7530 36.6 


0.1475 0.8708 -0.9029 18.1 
0.1845 1.1349 -0.9232 12.4 



4 i—i — . — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — r 

. WMAP 5 yr + SDSS BAO 2009 + 69 GRBs 

+ SHOES H + 397 SNe la ("Constitution" set) 




- 

. £ SNe la flux-averaged with dz = 0.03 at z^0.2 
cj> SNe la with no flux-averaging 



_j I i i i I i i i I i i i I i i i I i i i U 

0.2 0.4 0.6 0.8 1 

Z 

FIG. 5: Dark energy density function X(z) = px (z)/px(0) mea- 
sured from combining SN la data (the "Constitution" set of 397 
SNe la) with CMB, BAO, GRB data, and imposing the SHOES 
prior on Ho- The 68% (shaded) and 95% confidence level regions 
are shown. A flat universe is assumed. 



probability density distributions of the dark energy and 
cosmological parameters using [R(z*), l a {z*), z*, r s (z ( i)] 
(as we have done throughout this paper), and using 
[R(z*), l a (z*),u>t,]- Clearly, adding CMB constraints on 
r s {zd) tightens the constraints on (wq, w a ) and (wq, UJ0.5), 
as these are parameters from assuming a dark energy 
equation-of-state linear in the cosmic scale factor a, thus 
imply strong assumptions about early dark energy that 
are propagated to the low and intermediate redshifts [13] • 
We find that the constraints on 
[R(z*), l a {z*), 2*, fg^za)] summarize all useful infor- 
mation from CMB data that are relevant to dark 
energy and are geometric and independent of dark 



energy perturbations. Further adding constraints on 
Zd does not add new information, since Zd is well 
approximated by Eq. fTDj) , which only depends on Q m h 2 
and flbh 2 , and which in turn are already constrained by 
[R(z*),l a (z*), z*,r a (zd)]- 



IV. SUMMARY AND DISCUSSION 

It is likely that current sets of SNe la may be contam- 
inated by unknown systematic effects |43| . We have 
used the flux-averaging of SNe la (which was developed 
to reduce the weak lensing systematic effect [H], 5H, |46[ ) 
to help reduce the impact of unknown systematic effects. 
Since the SNe la in each redshift bin are flux-averaged 
and not used directly in the likelihood analysis, the sys- 
tematic bias that results from unknown systematic errors 
from individual SNe la can be minimized. For homo- 
geneous data without unknown systematic errors, flux- 
averaging should not change the results qualitatively at 
intermediate redshifts (where the lensing effect is negli- 
gible). Thus in addition to reducing lensing and lensing- 
like systematic effects, flux-averaging of SN la data pro- 
vide a useful test for the presence of unknown systematic 
effects at intermediate redshifts. 1 

We find that flux-averaging has a significantly smaller 
effect on the nearby+SDSS+ESSENCE+SNLS+HST 
data set of 288 SNe la [36 1, c ompared to the "Consti- 
tution" set of 397 SNe Ia[35[, see Figs ED and g] This 
same trend continues when the SN la data are combined 
with CMB, BAO, GRB data, and the SHOES measure- 
ment of Hq to measure the dark energy density func- 
tion, see Figs[3]|ni Furthermore, flux-averaging brings 



Note that flux-averaging corrects the bias in estimated distances 
due to weak lensing; this makes the statistical errors bigger when 
the number of SNe la flux-averaged in a bin is small. 
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FIG. 6: Dark energy density function X(z) = 

Px( z )/px (") measured from combining SN la data 
(nearby+SDSS+ESSENCE+SNLS+HST data set of 288 SNe 
la) with CMB, BAO, and GRB data, and imposing the SHOES 
prior on Hq. The 68% (shaded) and 95% confidence level 
regions are shown. A flat universe is assumed, (a) Assuming 
X(z > 1) = X[z = 1). (b) Same as (a), but with X[z = 1.5) 
added as a parameter, and assuming X(z > 1.5) = X (z = 1.5). 



FIG. 7: The measured cosmic expansion history H(z) correspond- 
ing to Fig[3] The error bars represent the 68% confidence level. 
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FIG. 8: The measured cosmic expansion history H(z) correspond- 
ing to Fig|4] The error bars represent the 68% confidence level. 



both data sets closer to the predictions of a flat uni- 
verse dominated by a cosmological constant (see FigsQJ 
[HI). The nearby+SDSS+ESSENCE+SNLS+HST data 
set of SNe la (together with CMB, BAO, GRB data, 
and the SHOES measurement of Ho) are consistent with 
a cosmological constant at 68% confidence level with 
or without flux-averaging of SNe la (see FigH]). With- 



out flux-averaging, the combined data using the "Con- 
stitution" set of 397 SNe la seem to indicate a devia- 
tion from a cosmological constant at ^95% confidence at 
^ z ^ 0-8, but is consistent with a cosmological con- 
stant at ~68% confidence when SNe la are flux-averaged 
(see FigE]). 
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WMAP5 [fl.l ,z.,r s (zj] + SDSS BAO 2009 + 69 GRBs 
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WMAP5 [fl,l ,z.,r s (zj] + SDSS BAO 2009 + 69 GRBs 
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FIG. 9: The 68% and 95% joint confidence level contours on con- 
straints on (wo, Wa) and (u>0, MJ0.5), using SN la data (the "Consti- 
tution" set of 397 SNe la) together with CMB, BAO, GRB data, 
and imposing the SHOES prior on Hq. 



FIG. 10: The 68% and 95% joint confidence level contours on 
constraints on (wo,w a ) and (wO, wo. 5), using SN la data (the 
ncarby+SDSS+ESSENCE+SNLS+HST set of 288 SNe la) to- 
gether with CMB, BAO, GRB data, and imposing the SHOES 
prior on Hq. 



It is interesting to note that for a dark energy den- 
sity function X(z) parametrized by {X(zi)} (^=0.25, 

0. 5, 0.75, 1.0) and assumed to be constant at z > 

1, the factor of improvement in the dark energy FoM 
is 5.2 when the nearby+SDSS+ESSENCE+SNLS+HST 
data set of 288 SNe la is used instead of the "Consti- 
tution" set of 397 SNe la (see Table VI). For a lin- 
ear dark energy equation of state, the improvement in 
FoM is a little less than ^50%. This indicates that 
current data already contain more information about 
dark energy than can be adequately represented by 
a linear dark energy equation of state. Using the 
nearby+SDSS+ESSENCE+SNLS+HST data set of 288 
SNe la gives more stringent constraints on dark energy 



than using the "Constitution" set of 397 SNe la, because 
the former data set includes 103 SNe la from SDSS at 
0.04 < z < 0.42 [36| which are not in the latter data set, 
and because the latter data set is less homogeneous and 
likely more prone to unknown systematic effects. 

We have measured the dark energy density function 
X(z) from data without imposing X(z) > 0. Since the 
probability density distributions for X(z) that we have 
obtained from the current data extend to X(z) < 0, it 
is not possible to convert the measurements of X(z) into 
measurements of wx(z) = px(z)/ Px( z ), or the deceler- 
ation parameter q(z) (which depends on wx(z)). X(z) 
is related to the dark energy equation of state wx{z) as 
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288 SNe la (flux-averag. dz=.03 at zi.2) + WMAP5 + SDSS BAO 2009 + 69 GRBs + SHOES H 



WMAP5 [R.l a ,z..r 1 (z <1 )] (solid); [R.l a ,uJ (dashed) 




-1 -0.5 -1.5 -1 0.25 0.3 

w o w o.6 n m 




0.05 0.65 0.7 0.75 0.8 0.022 0.024 



FIG. 11: Marginalized probability distributions of the dark 
energy and cosmological parameters using SN la data (the 
nearby+SDSS+ESSENCE+SNLS+HST set of 288 SNe la flux- 
averaged at z > 0.2) together with CMB, BAO, GRB data, and im- 
posing the SHOES prior on Ho, assuming a linear dark energy equa- 
tion of state parametrized by (wo, wo. 5) (see Eq.|20|). The solid 
and dashed lines indicate the use of the CMB data summarized by 
[R(z*),l a (z t ),z*,ra(Zd)], and [R(z,),l a (z,),uj b ] respectively. 

288 SNe la (inc. SDSS; dz=.03 flux-averag. at zit.2) + WMAP5 + SDSS BAO 2009 + 69 GRBs + SHOES H 
WMAP5 [R.l,.z..r,(z a )] (solid); [R.1..M,] (dashed) 
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FIG. 12: Same as Fig llll except the linear dark energy equation 
of state is parametrized by the usual (wo,w a ). 

follows 0: 

X{z) = MO) = p U dz i + z < / ■ (24) 

Hence parametrizing dark energy with wx{z) implicitly 



assumes that px{z) does not change sign in cosmic time 
(i.e., X(z) > 0). This precludes whole classes of dark 
energy models in which px(z) becomes negative in the 
future ("Big Crunch" models, see Refs. (l2l. Ifrij for an 
example) |57| . Thus the measurement of X(z), without 
assuming X{z) > 0, contains more information about 
dark energy than that of wx(z) or q(z). 

Our work complements Genovese et al. (2009) |65j . 
and Bogdanos & Nesseris (2009) |66[, both of these pa- 
pers assumed a flat universe, and analyzed SN la data 
only . Genovese et al. (2009) used model-testing (see e.g. 
|67j ]) to diffentiate assumptions about wx(z), and de- 
rived parametric and non-parametric estimators of wx{z) 
[65j . Bogdanos & Nesseris (2009) used genetic algorithms 
to derive model-independent constraints on the distance 
modulus from SNe la, which are then used to constrain 
wx(z) (66j . In this paper, we have derived model- 
independent constraints on X[z), instead of w x {z), for 
both a flat universe and a universe with arbitrary spatial 
curvature. We also use newer SN la data (which signifi- 
cantly tighten the dark energy constraints), and combine 
with CMB, BAO, GRB data, and the SHOES measure- 
ment of Hq. 

Finally, we note that the (wq, w a ) FoM for current data 
(the nearby+SDSS+ESSENCE+SNLS+HST data set of 
288 SNe la flux-averaged, with CMB, BAO, GRB data, 
and the SHOES prior on Ho) is 18.1, more than a fac- 
tor of two larger than that found by Wang (2008a) using 
data available in March 2008 (182 SNe la compiled by 
Riess et al. 2007 [H], [R(z»), l a {z*), fl b h 2 } from the five 
year WMAP observations, the SDSS measurement of the 
BAO scale by Eisenstein et al. 2005 [69], and assuming 
the HST prior of H = 72 ± 8 kms^Mpc" 1 ). Most of 
this gain results from the significantly improved distance 
measurements from SNe la (see Figsfl][2]). Assuming a 
flat universe, dark energy is detected at > 98% confidence 
level for z < 0.75 using the combined data with 288 SNe 
la from nearby+SDSS+ESSENCE+SNLS+HST, inde- 
pendent of the assumptions about X(z > 1) (see FigJB]). 

In order to make solid progress on probing dark en- 
ergy, it will be essential to launch aggressive observa- 
tional projects to obtain large and uniform sets of SNe 
la spanning the redshift range of < z < 2 [7(| • 
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